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Abstract. Inspired by the works of Rodnianski and Schlein [3T] 
and Wu [34l [35] , we derive a new nonlinear Schrodinger equation 
that describes a second-order correction to the usual tensor prod- 
uct (mean-field) approximation for the Hamiltonian evolution of 
a many-particle system in Bose-Einstein condensation. We show 
that our new equation, if it has solutions with appropriate smooth- 
ness and decay properties, implies a new Fock space estimate. 
We also show that for an interaction potential v{x) = ex{x)\x\~^ , 
where e is sufficiently small and x € even, our program can 
be easily implemented locally in time. We leave global in time 
issues, more singular potentials and sophisticated estimates for a 
subsequent part (part II) of this paper. 

1. Introduction 

An advance in physics in 1995 was the first experimental observation 
of atoms with integer spin (Bosons) occupying a macroscopic quantum 
state (condensate) in a dilute gas at very low temperatures [HH]. This 
phenomenon of Bose-Einstein condensation has been observed in many 
similar experiments since. These observations have rekindled interest 
in the quantum theory of large Boson systems. For recent reviews, see 
e.g. [231129]. 

A system of N interacting Bosons at zero temperature is described by 
a symmetric wave function satisfying the iV-body Schrodinger equation. 
For large N, this description is impractical. It is thus desirable to 
replace the many-body evolution by effective (in an appropriate sense) 
partial differential equations for wave functions in much lower space 
dimensions. This approach has led to "mean-field" approximations 
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in which the single particle wave function for the condensate satisfies 
nonlinear Schrodinger equations (in 3 + 1 dimensions). Under this 
approximation, the A^-body wave function is viewed simply as a tensor 
product of one-particle states. For early related works, see the papers 
by Gross [151 [I6], Pitaevskii |28l and Wu [Ml [35] . In particular, Wu [Ml 
[35] introduced a second-order approximation for the Boson many-body 
wave function in terms of the pair- excitation function, a suitable kernel 
that describes the scattering of atom pairs from the condensate to other 
states. Wu's formulation forms a nontrivial extension of works by Lee, 
Huang and Yang [21j for the periodic Boson system. Approximations 
carried out for pair excitations [2T1 [Ml [3S] make use of quantized fields 
in the Fock space. (The Fock space formalism and Wu's formulation 
are reviewed in sections ll.ll and [1.31 respectively.) 

Connecting mean-field approaches to the actual many-particle Hamil- 
tonian evolution raises fundamental questions. One question is the 
rigorous derivation and interpretation of the mean field limit. Elgart, 
Erdos, Schlein and Yau [HI [3 [SI [3 [101 E] showed rigorously how mean- 
field limits for Bosons can be extracted in the limit — > cxd by using 
Bogoliubov-Born-Green-Kirkwood-Yvon (BBGKY) hierarchies for re- 
duced density matrices. Another issue concerns the convergence of the 
microscopic evolution towards the mean field dynamics. Recently, Rod- 
nianski and Schlein [31] provided estimates for the rate of convergence 
in the case with Hartree dynamics by invoking the formalism of Fock 
space. 

In this paper, inspired by the works of Rodnianski and Schlein [31] 
and Wu [Ml [M] , we derive a new nonlinear Schrodinger equation de- 
scribing an improved approximation for the evolution of the Boson 
system. This approximation offers a second-order correction to the 
usual tensor product (mean field limit) for the many-body wave func- 
tion. Our equation yields a corresponding new estimate in Fock space, 
which complements nicely the previous estimate [31]. 

The static version of the many-body problem is not studied here. The 
energy spectrum was addressed by Dyson [5] and by Lee, Huang and 
Yang [21]. A mathematical proof of the Bose- Einstein condensation 
for the time-independent case was provided recently by Lieb, Seiringer, 
Solovej and Yngvanson [221 [251 [2il [25] . 

1.1. Fock space formalism. Next, we review the Fock space over 
L^(M'^), following Rodnianski and Schlein [31] • The elements of are 
vectors of the form '0 = {ipo, %Iji[xi) ,'il)2{xi, X2) ^ ■ ■ ■) , where ipQ & and 
ipn G L^(]R^'^) are symmetric in xi, . . . , x„. The Hilbert space structure 
of is given by (0, ^l^) = Y.n I (f^n'ipndx. 



SECOND ORDER CORRECTIONS , I 



3 



For / G L^(M^) the (unbounded, closed, densely defined) creation 
operator a*{f) : T ^ T and annihilation operator a(/) : T ^ T are 
defined by 



1 " 

.7=1 



(a*(/)^n-l) [Xx,X2, 
{a(f)tlJn+l) {Xi,X2, ■■■ ,Xn) = Vn+1 j ^(n+l){x,Xi, ■ ■ ■ ,Xn)fix) dx . 

The operator valued distributions a* and defined by 

a* if) = / fix)al dx , 



a{f) = J f{x) dx . 

These distributions satisfy the canonical commutation relations 

[a^,a*] =6{x-y) , (1) 
[a^, ay] = [<, a*] = . 

Let A?" be a fixed integer (the total number of particles), and v{x) be 
an even potential. Consider the Fock space Hamiltonian : T ^ T 
defined by 

Hn^ J alAa^dx+^ j v{x - y)ala*ya^ay dxdy (2) 

This Hn is a diagonal operator which acts on each in correspondence 
to the Hamiltonian 

n ^ 
HN,n = J2 + ^ Y^ ^i^i - ^i) ■ 

In the particular case n — N, this is the mean field Hamiltonian. 
Except for the introduction, this paper deals only with the Fock space 
Hamiltonian. The reader is alerted that "PDE" Hamiltonians such as 
HN,n will always have two subscripts. The sign of v will not play a 
role in our analysis. However, the reader is alerted that due to our 
sign convention, i> < is the "good" sign. The time evolution in the 
coordinate space for Bose- Einstein condensation deals with the function 

e^*-^"'"V'o (3) 



4 



M. GRILLAKIS, M. MACHEDON, AND D. MARGETIS 



for tensor product initial data, i.e., if 

where ||0o||l2(k3) = 1. This approach has been highly successful, even 
for very singular potentials, in the work of Elgart, Erdos, Schlein and 
Yau O [3 El 13 ini HI]- In this context, the convergence of evolution 
to the appropriate mean field limit (tensor product) as — > oo is 
established at the level of marginal density matrices 'j^^^ in the trace 
norm topology. The density matrices are defined as 

1.2. Coherent states. There are alternative approaches, due to Hepp 
|17j . Ginibre and Velo [13], and, most recently, Rodnianski and Schlein 
|31j which can treat Coulomb potentials v. These approaches rely on 
studying the Fock space evolution e^^^^xf^Q where the initial data tpQ is 
a coherent state, 

■00 = (co, ci0o(a;i), C20o(a;i)0o(a;2), ■ ■ ■ ) ; 

see dl]) below. The evolution ^ can then be extracted as a "Fourier 
coefficient" from the Fock space evolution; see [3^. Under the as- 
sumption that f is a Coulomb potential, this approach leads to strong 
L^-convergence, still at the level of the density matrices '^\^\ as we 
will briefly explain below. 

To clarify the issues involved, let us consider the one-particle wave 
function (j){t, x) (to be determined later as the solution of a Hartree 
equation), satisfying the initial condition 0(0, a;) = (t)Q{x). Deflne the 
skew-Hermitian unbounded operator 

A(0) = a(0)-a*(0) 

and the vacuum state VL = (1, 0, 0, ■ ■ ■ ) G JF. Accordingly, consider the 
operator 

ty(0) = e-v^^W , 

which is the Weyl operator used by Rodnianski and Schlein [21]. The 
coherent state for the initial data 00 is 

= e-^MV2 1^1^ ... ^ 1^^^ . . . ... V (4) 
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Hence, the top candidate approximation for e'^^^Vo reads 

^tensor(i) = e-^^(^(*'-))fi . (5) 

Rodnianski and Schlein [31] showed that this approximation works (un- 
der suitable assumptions on v), in the sense that 



iTr 



pCt 

= 0(— ) iV^oo; 

the symbol Tr here stands for the trace norm in x G M"' and y G 
M^. The first term in the last relation, including -j^, is essentially the 

density matrix 'y[^\t,x,y). For the precise statement of the problem 
and details of the proof, see Theorem 3.1 of Rodnianski and Schlein 

m- 

Our goal here is to find an explicit approximation for the evolution 
in the Fock space. For this purpose, we adopt an idea germane to 
Wu's second-order approximation for the A^-body wave function in Fock 
space [23 [35]. 

1.3. Wu's approach. We first comment on the case with periodic 
boundary conditions, when the condensate is the zero-momentum state. 
For this setting, Lee, Huang and Yang [21] studied systematically the 
scattering of atoms from the condensate to states of opposite momenta. 
By diagonalizing an approximation for the Hamiltonian in Fock space, 
these authors derived a formula for the A^-particle wave function that 
deviates from the usual tensor product, as it expresses excitation of 
particles from zero monentum to pairs of opposite momenta. 

For non-periodic settings, Wu [3ll L35J invokes the splitting = 
ao{t)4>{t, x) + ax,i{t) where oq corresponds to the condensate, [oq, ag] = 
1, and ttxA corresponds to states orthogonal to the condensate, [oq, = 
= [ao,a*. ]^]. Wu applies the following ansatz for the A^-body wave 
function in Fock space: 

Ar(t)e^[^«V^(t) , (6) 

where il>%{t) describes the tensor product, A/'(t) is a normalization fac- 
tor, and P[-ft'o] is an operator that averages out in space the excitation 
of particles from the condensate to other states with the effective 
kernel (pair excitation function) Kq. An explicit formula for P[i^^o] is 



V[Ko] = [2No{t)]-' / al^.al, Ko{t, x, y) a,{tY , (7) 
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where A^o is the expectation value of particle number at the condensate. 
This Kq is not a-priori known (in contrast to the case of the classical 
Boltzmann gas) but is determined by means consistent with the many- 
body dynamics. In the periodic case, ([6]) reduces to the many-body 
wave function of Lee, Huang and Yang (21j . 

Wu derives a coupled system of dispersive hyperbolic partial differen- 
tial equations for (0, Kq) via an approximation for the A^-body Hamil- 
tonian that is consistent with ansatz A feature of this system is 
the spatially nonlocal couplings induced by Kq. Observable quantities 
such that the depletion of the condensate can be computed directly 
from solutions of this PDE system. This system has been solved only 
in a limited number of cases [351 ISSl |2Z] ■ 

1.4. Scope and outline. Our objective in this work is to find an 
explicit approximation for the evolution 

in the Fock space norm, where i/'o is the coherent state (jl]). This would 
imply an approximation for the evolution 

in L^(M^^) as N ^ oo. To the best of our knowledge, no such ap- 
proximation is available in the mathematics or physics literature. In 
particular, the tensor product type approximation ^ for satisfying 
a Hartree equation, as in [3T], is not known to be such a Fock space 
approximation (nor do we expect it to be). 

To accomplish our goal, we propose to modify ([5]) in two ways. One 
minor correction is the multiplication by an oscillatory term. A second 
correction is a composition with a second-order "Weyl operator" . Both 
corrections are inspired by the work of Wu [511 ES]; see also 1^ . 
However, our set-up and derived equation is essentially different from 
these works. 

We proceed to describe the second order correction. Let k{t,x,y) = 
k(t, y, x) be a function (or kernel) to be determined later, with fc(0, x, y) = 
0. The minimum regularity expected of /c is A; G L^{dxdy) for a.e. t. 

We define the operator 

i? = -y (^k{t, x,y)aj.ay — k(t, x,y)ala*) dx dy . (8) 

Notice that B is skew-Hermitian, i.e., iB is self-adjoint. The operator 
could be defined by the spectral theorem; see [30]. However, we 
prefer the more direct approach of defining it first on the dense subset 
of vectors with finitely many non-zero components, where it can be 
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defined by a convergent Taylor series if ||fc||L2(f;2;^y) is sufficiently small. 
Indeed, B restricted to the subspace of vectors with all entries past the 
first identically zero has norm < CA^||A;||i2. Then is extended to 
as a unitary operator. 

Now we have described all ingredients needed to state our results 
and derivations. The remainder of the paper is organized as follows. 
In section [2] we state our main result and outline its proof. In section [3] 
we study implications of the Hartree equation satisfied by the one- 
particle wave function (j){t,x). In section H] we develop bookkeeping 
tools of Lie algebra for computing requisite operators containing B. In 
section Owe study the evolution equation for a matrix K that involves 
the kernel k. In section |6] we develop an argument for the existence of 
solution to the equation for the kernel k. In section [7]we find conditions 
under which terms involved in the error term e^Ve~^ are bounded. In 
section [8] we study similarly the error term e^[A, V]e~^. In section |9] 
we show that we can control traces needed in derivations. 



2. Statement of main result and outline of proof 

In this section we state our strategy for general potentials satisfying 
certain properties. Later in the paper we show that all assumptions of 
the related theorem are satisfied locally in time for v{x) = x{^)'^\y 
sufficiently small, and x ^ C^: even. 

Theorem 2.1. Suppose that v is an even potential. Let (p be a smooth 
solution of the Hartree equation 

«^ + A0+(t;*|0|')0 = O (9) 

with initial conditions (pQ, and assume the three conditions listed below: 

(1) Assume that we have k{t,x,y) G L'^{dxdy) for a.e. t, where k 
is symmetric, and solves 

{iut + ug^ + gu- {I + p)m) = {ipt + [g,p] + um){l + p^^u , (10) 

where all products in ( JTOl) are interpreted as spatial compositions 
of kernels, "1 " is the identity operator, and 

u(t,x,y) := sh.{k) := k + —kkk + . . . , (11) 

6{x-y)+ p{t, X, y) := ch{k) ■= 6{x - y) + ^kk + ... , 

git, X, y) := -AJ{x - y) - v{x - y)(f){t, x)0(t, y) - (v * |0n(t, x)6{x - y) , 
m(t, X, y) := v{x - y)4>{t, x^it, y) . 
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(2) ^4/50, assume that the functions 

fit) := \\e''[A,V]e-''n\\:r 

and 

git) := We^'Ve-^ny 

are locally integrable (V is defined in ([2])j. 

(3) Finally, assume that J dit, x, x) dx is locally integrable in time, 
where 

dit, X, y) = {ishik)t + shik)g'^ + gshik)) shjk) 
-itchik)t + [g,chik)])diik) 



— sh(/c)mch(/c) — ch(/c)msh(A;) . 
Then, there exist real functions Xo, Xi such that 

^ lo fis)ds jQgis)ds 



Recall that we defined (see sectionU^ 

■00 = e~^^^^^^l an arbitrary coherent state (initial data) , 
Ait) = ai^it,-))-a*i<Pit,-)) , 

Bit) ~ 2 J {^{^■'^■>y)^xa.y ~ ^i^->^->y)^*x^*y) dxdy . 
A few remarks on Theorem 12.11 are in order. 
Remark 2.2. Written explicitly, the left-hand side of (fTOl) equals 

iut + ug'^ + gu- il+ p)m = - - A^^ uit, x, y) 



— (f)it,x) J vix — z)(j)it, z)uit, z,y) dz — (f)it,y) J uit,x, z)viz — y)(j)it, z) dz 

-iv* x)uit, X, y)-iv* {(pD^it, y)uit, x, y) 

- vix - y)4>it,x)'4>it,y) 

-'(P{'t,y) / il +p)it,x,z)viz - ?/)0(t,^) dz . 
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The main term in the right-hand side equals 

_ d 

iPt + [g, p]+um = i-g^p{t, X, y) + (-A^ + Aj^) p(t, x, y) 



— (f){t,x) J v{x — z)(j){t, z)p{t, z,y) dz 

+ <P{t,y) j pit,x,z)v{z - y)'^{t,z) dz 

-{v* |0n(t, x)p{t, X, y) + {v* \4>\ f{t, y)p{t, x, y) 

+ u{t, X, z)v{z — y)(f){t, z)(l){t, x) dz . 



Remark 2.3. The algebra, as well as the local analysis presented in this 
paper do not depend on the sign of v. However, the global in time 
analysis of our equations would require v to be non-positive. 

Remark 2.4. Our techniques would allow us to consider more general 
initial data of the form i/?q = e'^^^^^e^^'^^^Vl. For convenience, we 
only consider the case of tensor products (-B(O) = 0) in this paper. 

Proof. Since and are unitary, the left-hand side of ( fT2l) equals 

^\QiIoi^Xois)+Xiis))ds^Bit)^VNA{t)^itHN^-VNA{0)^ _ ^||^ _ 

Define 

$(t) = eBit)^VNAit)^itH^-^A{0)^ 

In Corollary 15.21 of section Owe show that our equations for (p, k insure 
that 

1 d 

--vI/ = Lv[/ , 

I ot 

where L = L — Nxo — Xi fo^^ some L: Hermitian, i.e. L = L*, where 
L commutes with functions of time, Xoy Xi are real functions of time, 
and, most importantly (see corollary 15.21 of section [5] and the remark 
following it), 

\\Ln\\^ < N-^/^\\e''[A, Vje'^'nil^ + N-'\\e''Ve-''n\\^ . (13) 
We apply energy estimates to 

(1^ _ 1^ ^e^foi^xo{s)+xi(s))ds^ -il)=Ln. 
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Explicitly, 



< 2 {N-^/^\\e^[A,V]e-^n\\^ + N-'\\e^Ve~^n\\^) || (e^/o(^xoW+xi)rf^^ _ 
Thus 

^||(e^/o(^xoW+xi)<is^ _^)|| < N-'/^\\e^[A,V]e~^Q\\:F + N-^\\e^Ve-^n 

and (fT2l) holds. This concludes the proof. □ 
□ 

3. The Hartree equation 

In this section we see how far we can go by using only the Hartree 
equation for the one-particle wave function (p. 

Lemma 3.1. The following commutation relations hold (where the t 
dependence is suppressed, A denotes A{(j)) and V is defined by formula 



[A, V] = j vix - y) {(l){y)ala^ay + (p{y)ala*ya,^) dxdy 
[A, [A, V]] (14) 
v{x - y) [(j){y)(j){x)axay + (f){y)(j){x)ala*y + 2(j){y)(f){x)alay) dxdy 

+ 2 J [v * 10^1) {x)al.ax dx 

A [A,[AV]] 

= 6 y [v * 10^1) (x) (0(x)a* + 0(x)aa.) dx 

A A[A,[AV]] 
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Proof. This is an elementary calculation and is left to the interested 
reader. □ 

Now, we consider \E'i(t) = e^^'^*)e**^e~^^*^°''fi for which we have 
the basic calculation in the spirit of Hepp [T7] , Ginibre- Velo [13] , and 
Rodnianski-Schlein pT]; see equation (3.7) in [31] . 



Proposition 3.2. If (p satisfies the Hartree equation 

+ A<f) + {v * \<f)\^)<P = 

while 



then \l'i(t) satisfies 



~idi 



+N-^'^[A, r] + iV-V-y j v{x ~ y)\ct>{t,x)\''\(t>{t,y)\Hxdy^^!^{t) 



Proof. Recall the formulas 



e^iJe-^ = H+[C,H] + ^^ [C, [C, H]] + 



and 



Applying these relations to C = v NA we get 

1 d 
I ot 



(15) 



where 



i \dt J 
i (^N^/^A + ^[A,A]^+H + iVV2[A, Ho] 



+ N-'"[A,V] + ^[A,[A,H, 



+ 



N 
4! 



A, 
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Eliminating the terms with a weight of V^, or setting 

, (16) 



U+[A,Ho] + ^U, [A, [A,V]\ 



is exactly equivalent to the Hartree equation ([9]). By taking an addi- 
tional bracket with A in f|T6|) . we have 



^[A,A]+[A, [A,Ho]]+^^ 



A, 



A, [A, [A,V]] 



0, 



and thus simplify (ITSl) to 
1 d 



^g^Mt)= \H, + -[A,[A,V]] 



+N-^/^[A, V] + N-^V -N 



4! 



A, 



A, [A, [A, V]] 



□ 



This concludes the proof. 

□ 

The first two terms on the right-hand side are the main ones. The 
next two terms are O and O (^j^). The last term equals 

N f 

-y / v{x - y)\(l){t,x)\^\(f){t,y)\'^dxdy:= -Nxo ■ 

Notice that ||Li(f2)|| is not small because of the presence of a*a* in 
[a, [A, V]\. In order to eliminate these terms, we introduce B (see (|H])) 
and take 



Accordingly, we compute 



where 



1^ / T I 

I ot 



and 



1 f d 



i \dt 



Lq + Ar-i/2e^[A, y]e-^ + N-'e^'Ve-'' - Nxo 



contains all quadratics in the operators a, a*. 



(17) 
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Equation ( ITOi) for k turns out to be equivalent to the requirement 
that L has no terms of the form a* a* . Terms of the form aa* will 
occur, and will be converted to a* a at the expense of Xi- 

In other words, we require that Lq have no terms of the form a*a* . 
For a similar argument (but for a different set-up), see Wu [35] . 



4. The Lie algebra of "symplectic matrices" 

In this section we describe the bookkeeping tools needed to compute 
Lq of f|T7j) in closed form. The results of this section are essentially 
standard, but they are included here for the sake of completeness. 

We start with the remark that 

Wi) + a*{g,),a{h) + a*{g2)] = J fig2 - f2gi (18) 

where 

j_f -5ix - y)\ 
- [5ix -y) )■ 

This observation explains why we have to invoke symplectic linear al- 
gebra. We thus consider the infinite-dimensional Lie algebra sp of 
"matrices" of the form 



for symmetric kernels k = k{t,x,y) and / = l{t,x,y), and arbitrary 
kernel d{t,x,y). (The dependence on t will be suppressed when not 
needed.) This situation is analogous to the Lie algebra of the finite- 
dimensional complex symplectic group, with x, y playing the role of i 
and j. We also consider the Lie algebra Quad of quadratics of the form 

Q(d,k,l):J-(a, al)(l (^^J) (19) 



a^a* + a*ax 1 
d{x,y) ^-^ — - — dxdy + - j k{x,y)axay dx dy 



^ j Kx,y)(i*xO*y dxdy 



{k, I and d as before). Furthermore, we agree to identify operators 
which differ (formally) by a scalar operator. Thus, J d{x,y)axa* is 
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considered equivalent to / d{x,y)a*ax- We recall the following result 
related to the metaplectic representation (see, e.g. [T2]). 

Theorem 4.1. Let S = S{d, k,l), Q = Q{d, k, I) related as above. Let 
f , g be functions (or distributions). Denote 

:= / ifix)a^ + g{x)al)dx . 



k9. 

We have the following commutation relation: 

[Q,K,<)(^)] = K,<)5(^) (20) 

where products are interpreted as compositions. We also have 

e^K,<)(^)e-« = («.,«:)e"(^) , (21) 

provided thate^ makes sense as a unitary operator (Q: skew-Hermitian) . 

Proof. The commutation relation (!20l) can be easily checked directly, 
but we point out that it follows from (fTSl) . In fact, using (JT8l) . for any 
rank one quadratic we have 

[(a(/i) + a*{g,)) (aif^) + a\g^)) , a(/) + a\g)] 

7' 



Thus, for any R we have 

Now specialize to i? = \SJ, S G sp, and use 5*^ = JSJ to complete 
the proof. 

The second part, equation fl2T|) . follows from the identity 



[(a. al)R(^^)j,a{f) + a*{g)] = -{a^ al){R + R^)j 



eQce-« = c + [g,c] + ^[g,[g,c]] +... , 

or, in the language of adjoint representations, Ad(e'^)(C) = e^*^*^'^^(C), 
which is applied to C = a(/) + a*{g). □ 

□ 

A closely related result is provided by the following theorem. 
Theorem 4.2. (1) The linear map I : sp ^ Quad defined by 



S{d,k,l) Q{d, k,l) 
is a Lie algebra isomorphism. 
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(2) Moreover, if S = S{t), Q = Q{t) andl{S{t)) = Q{t) is skew- 
Hermitian, so that is well defined, we have 



(3) Also, if R E sp, we have 

I {e^Re-^) = e^I{R)e~'^ . (23) 

Remark 4.3. In the finite-dimensional case, this is (closely related to) 
the "infinitesimal metaplectic representation"; see p. 186 in [12] . In 
the infinite dimensional case, we must be careful, as some of our oper- 
ators are not of trace class. For instance, J a^a* does not make sense. 

Proof. First, we point out that implies (12^ . at least in the case 
where R is the "rank one" matrix 



Notice that fl2Tl) can also be written as 



In conclusion, we find 
[ax a*) R 



ay 



a 



X 



e« (a. al) {^^ e~^e^ [h i) J 
(a. al) {h i) Je'^^ 



[a^ a*) e^Re' 



-s -a; 



ay 



since S"^ = JSJ ii S ^ sp, and Je-^^-^ = e^^J. 

We now give a direct proof that fll9p preserves Lie brackets. Denote 
the quadratic building blocks by Qxy = a^ay, Qly = a^a*, N^y = 
I (^a^a* + a* ax) ■ One can verify the following commutation relations. 
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which will be also needed below: 



[Qxy, Qlw] =S{x - z)Ny^ + 5{x- w)Ny^ + 6{y 

[Qxy, N^y,] =6{x - w)Qyz + S{y - w)Q^^ , 
[Nxy, QL] =S{x- z)Q*y^ + 6{x- w)Qy-, , 
[N^y, N:,J\ =6{x - w)N^y - 6{y - z)N^^ . 
Using ^2M we compute 

1 
2 



k{x^y)axaydxdy, / l{x,y)al.a*dxdy 



which corresponds to the relation 

k 





/ 



kl 





-Ik 



(24) 
(25) 
(26) 
(27) 



{kl){x,y)Nxy dxdy 



The other three cases are similar. 

To prove fl22l) . expand both the left-hand side and the right-hand 
side as 



= Q + l[Q,Q] + --- 




The proof of 0231) is along the same lines. □ 
□ 

Remark 4.4. Note on rigor: All the Lie algebra results that we have 
used are standard in the finite-dimensional case. In our applications, S 
will be K where is a matrix of the form (|29l) . see below, and Q will 
be i? = 1{K). The unbounded operator B is skew-Hermitian and e^ip 
is defined by a convergent Taylor series iiip ^ ^ has only finitely many 
non-zero components, provided ||fc(t, •, •)llL2(dxdy) is small . We then 
extend to all as a unitary operator. The norm \\k{t, ■, ■)\\L'^(dxdy) 
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iterates under compositions; thus, the kernel is well defined by its 
convergent Taylor expansion. In the expression 

gSpg-B = p + p] + _ _ _ (28) 

for P, a first- or second-order polynomial in a, a*, we point out that the 
right-hand side stays a polynomial of the same degree, and converges 
when applied to a Fock space vector with finitely many non-zero com- 
ponents. For our application, we need to know ( |28l) is true when applied 
to VL. The same comment applies to the series 

5. Equation for kernel k 
Now apply the isomorphism of the previous section to the operator 

B = I{K) 

for 



This agrees to the letter with the isomorphism (]T9|) . The next two 
isomorphisms, (15U]) and ([SI]), require special treatment because aa* 
terms mirroring the a* a terms are missing in ([2]), (I14p . However, the 
discrepancy only happens on the diagonal. Once the relevant terms are 
commuted with B, they fit the pattern exactly. It isn't quite true that 



Ho=I 



-{A6){x-y) 

{M){x-y) 



-(("o" 1)) <») 

since, strictly speaking, 

iA6){x-y) \\ r a^Aa^^a^Aa^^^ 



{A6){x-y) 

is undefined. However, one can compute directly that [Axax,a*] = 
{A6){x — y) . Using that, we compute 

[S, Ho] = ^ j ((A^. + Ay)k{x, y)axay + (A^; + Ay)k{x, y)ala*y) dx dy . 

This commutator is in agreement with (l29l) . (!30|) . and the result can 
be represented in accordance with (fT9|l . namely 
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We also have 



k 
k 



iA6){x-y) 

{A6){x-y) 



X e 



-{A6){x-y) 

{A6){x-y) 



iA5){x-y) 

{A5){x-y) 



since c^Hqc ^ — Hq = [B, Hq] + | [B, Hq]^ + ■ ■ ■ The same comment 
apphes to the diagonal part of 



2 [A [A,V]] 



I 



12 



(31) 



-t;i20i02 - {V * 5i2 _ f 120102 

-Vi2(f)l4>2 ^"120102 + {V* 101 

where fi20i02 is an abbreviation for the product v{x — y)(f){x)(f){y), etc. 
Formula fl3T|) isn't quite true either, but becomes true after commuting 
with B. 

To apply our isomorphism, we quarantine the "bad" terms in (I5U1) 
and the diagonal part of ( l3Ti) . Define 



G 



9 




and 



M 



m 
— m 



where 



g = -A5i2 - ^^120102 - {v* 
m = f 120102 , 



12 



i/o + ^[A [A,V]] =Hg + I{M) 



and split 



where 

Hq = Hq+ J v{x — y)(f){y)(f){x)al.ay dxdy 

+ j [y * 10^1) {x)al.ax dx . 

By the above discussion we have 

[B,Hg]=I{[K,G]) and 
[e^,i/G]e-^=J([e^,G]e-^) . 



(32) 
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Write 

B 



I \dt 

+Hg + [e^, Hcje-"" + e«J(M)e-^ 

= HG+IiMi + M2 + M3) . (33) 
Notice that if K is given by ( l29l) . then 



e 

where 



^ ^ /ch(fc) sh(fc) \ 
Vsh(A;) ch(A:)>' 



ch(A;) = I + ^kk + ^kkkk + ... , (34) 

and similarly for sh(A;). Products are interpreted, of course, as compo- 
sitions of operators. 
We compute 



1 fdi{k)t sh(fc) A / ch(fc) -sh(fc )^ 
^ z Vsh(A;), ch(A;)J V-sh(A;) ch(A;) ^ 

1 / ch(A;)tch(A;) - sh{k)tsh(k) -ch{k)tsh{k) + sh(A;)tch(A;)\ 



^ /^[ch(A;),^] -sh{k)g'^ - gsh{k)\ 



and 



^2 = [e^,G]e-^ = 

[ ch, g] ch + (sh^f^ + (yfsh)sh — [ ch, g]sh — {shg'^ + 5fsh)ch\ 
* * y 

where sh is an abbreviation for sh(A;), etc, and 

-shm ch — chmsh shmsh + chmch \ 



= e^Me 
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Now define 

M = Ml +M2 + M3 . 
We have proved the following theorem. 

Theorem 5.1. Recall the isomorphism (fT9|) of Theorem \4.S\ 

(1) If Lq is given by (fT7|) . then 

Lq=Hq + j v{x - y)(p{y)(f){x)alay dx dy (35) 

+ J {v*\(f)^\) {x)ala^ dx + X{M) . 

(2) The coefficient of axdy in I (Ai) is —M.12 or 

{is\i{k)t + sh(fc)/ + ^sh(A;))ch(A;) - {ic\i{k)t - [ch(A;), ^])sh(A;) 



— sh(/c)msh(/c) — ch(/c)mch(fc) . 

(3) The coefficient of a*a* equals minus the complex conjugate of 
the coefficient of a^ay. 

(4) The coefficient of M-u, or 



d(t, X, y) = (^ish{k)t + sh{k)g^ + gsh{k)) sh{k) 
-{tch{k)t + [g,ch{k)])ch{k) 



—sh.{k)mch(k) — ch.{k)msh.{k) . (36) 

Corollary 5.2. If (p and k satisfy IQ and f fTOl) of theorem (12. ip . then 
the coefficients of a,j.ay and al.a* drop out and Lq becomes 

Lq =Hq + J v{x — y)(l){t^y)(f){t, x)al.ay dx dy + J [v * {(f)^]) {x)al.ax dx 

d{t,x,y) dxdy , 

where d is given by ( l36l) and the full operator reads 



L=Ho+ / v{x — y)(j){y)(j){t^ x)alaydxdy + / (t> * (x)a*aa;(ix 



j d{t, X, y)a*yaxdx + N-^''^e'^[A, V]e-^ + N~^e^Ve-^ - Nxo - Xi 



:= L - Nxo - Xi , 

and 



Xo = ^ / v{x -y)\(l){t,x)\'^\(f){t,y)\''^dxdy 
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= — - J d{t,x,x)dx . 

Remark 5.3. Notice that 

In = {N~'/^e^[A, V]e-^ + N-'e^Ve-^) , 

and therefore we can derive the bound 

\\Ln\\ < N-'^''\\e^[A,V]e-^n\\ + N-'\\e^Ve-^n\\ . 

Also, L is (formally) self-adjoint by construction. The kernel d{t,x,y), 
being the sum of the (1,1) entry of the self-adjoint matrices (j^e^) e~^, 
[e^, G]e~^ = e^Ge~^ —G and the visibly self-adjoint term — sh(fc)mch(fc)- 
ch(/c)msh(/c), is self-adjoint; thus, it has a real trace. Hence, L is also 
self-adjoint. 

In the remainder of this paper, we check that the hypotheses of our 
main theorem are satisfied, locally in time, for the potential v{x) = 

6. Solutions to equation ffTOl) 

Theorem 6.1. Let eo be sufficiently small and assume that v{x) = j^, 
orv{x) = x{^)j^\ forx ^ C^(]R^) . Assume thatcj) is a smooth solution 
to the Hartree equation (fT6ll . ||0||L2(rf^) = 1. Then there exists k G 
L°°([0, l])L'^{dxdy) solving ffTOj) with initial conditions k{0, x,y) = for 
< t < 1. The solution k satisfies the following additional properties. 

(1) 

i— — — Ay j /c||L°°[o,i]i2(dxdi/) < C . 



(2) 



(3) 



dt 

d 



— Ax — sh(A;)||ioo[o,i]L2(^^^y) < C 



i-Q^ — + Ay 1 p\\L°°[0,l]L^{dxdy) < C . 

(4) The kernel k agrees on [0, 1] with a kernel k for which 

\\k\\^i^..<C; 

see (138|) for the definition of the space X'^'^ and, of course, ^+ 
denotes a fixed number slightly bigger than \. 



S = i— - - Ay 



T = l—-A.,+ Ay 
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Proof. We first establish some notation. Let S denote the Schrodinger 
operator 

d_ 

'dt 

and let T be the transport operator 

d_ 
dt 

Let e : L'^{dxdy) L?{dxdy) denote schematically any linear operator 
of operator norm < Ceg, where C is a "universal constant". In practice, 
e will be (composition with) a kernel of the type </)(t, x)(j){t, y)v{x — y), 
or multiplication by f * Also, recall the inhomogeneous term 

m{t, X, y) = v{x - y)4>{t, x)4>{t, y) . 

Then, equation (ITU]) , written explicitly, becomes 

Sk = m + S{k - u) + e{u) + e{p) + [Tp + e{p) + e{u)){l + p)-\ . (37) 



Note that ch(fc)^ — sh.{k)sh.{k) = 1; thus, 1 + p = ch.{k) > 1 as an 
operator and (1 +p)~^ is bounded from to L^. We plan to iterate 
in the norm N{k) = || A;||L°=[o,i]L2{dxdy) + \\Sk\\L^io^i]L2(^axdy)- Notice that 

||'^||L2(da;d?/) < Ceo. 

Now solve 

Sko = m 

with initial conditions A;o(0, ■, ■) = 0, where N^ko) < Ceo. Define uq, po 
corresponding to fco- 

For the next iterate, solve 

Ski = m + S{ko-uo) + e{uo) + e{po) + (Tpo + e{po)+e{uo)){l+Po)~^uo ; 

the non-linear terms satisfy 

II ^(no — A;o)||Loo[o^i]i2(^^^y) = 

11^ ({Sko)koko - ko{Sko)ko + kohSko^ H \\L°-io,i]L\dxdy) 

= OiN{kof) . 
Also, recalling that po = ch(A;o) — 1, we have 

\\T{po)\\L°°[0,l]L^{dxdy) = II - ({Sko)ko — ko{Sko)^ + ■ ■ ■ \\L°°[0,l]L2{dxdy) 

= o{N{kor) . 

Thus, N{ki) < Ceo + Cel. Continuing this way, we obtain a fixed 
point solution in this space which satisfies the first three requirements 
of theorem 16.11 
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In fact, we can apply the same argument to j D'^k, since (^J^j V^m G 
L°°[0, l]L'^{dxdy) for < a < However, we cannot repeat the argu- 
ment for D^/'^k. 

We would like to have \\SD^'^k\\ 

L°°[o,i]L2(dzdj/) finite. Unfortunately, 
this misses "logarithmically" because of the singularity of v. 

Fortunately, we can use the well-known X'^'^ spaces (see |2l UHl [20] ) 
to show that || |>S'|*D"^/^M||i2(^j)2.2((ixd?;) is finite locally in time for (all) 
1 > s > \. This assertion will be sufficient for our purposes. Recall 
the definition of X^'^\ 

\m{\r - l^ri + lYu\Widrdi) ■■= ll^llxM . (38) 

Going back to (1371) . we write 

S{k) =m + F 

where we define the expression 

F{k) := S{k — u) — e{u) + pm + (T{p) + e{p) + um) (1 + p)~^u . 

The idea is to localize in time on the right-hand side: 

sCk)=xit){m + F) , 

where x^C^{^),X=^ on [0, 1]. Then, k = k on [0, 1]. 

As we already pointed out, we can estimate \\S (^) D"-k\\L2[o,i]L2{dxdy) < 
C for < a < i. We can further localize k in time to insure that 
these relations hold globally in time. By using the triangle inequality 
k^l'CPl + kl^ |^P)We immediately conclude that 

□ 

□ 



7. Error term e^Ve~^ 

The goal of this section is to list explicitly all terms in e^Ve~^ and 
to find conditions under which these terms are bounded. Recall that 
V is defined hj V = J v{xo - |/o)Q^oj/o^^oyo ^Xq dyo. For simplicity, 
shb(A;) denotes either sh(/c) or sh(A;), and chb(A;) denotes either ch(A;) 
or ch(A;). 

Let Xq ^ Uq] we obtain 



O* = J 
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According to the isomorphism ( |T9i) . we have 



-26{x - xo)6{y - yo) 

where the operator 

^ //ch(A;) sh(A;)\ / 0\ / ch(A;) -sh(fc)\\ 

~ Vvsh(^) ch(A;)y' V-25(a;-a;o)5(y-yo) Oyl V-sh(A;) ch{k) J J 
is a hnear combination of the terms 

chb(.)(.,.o)chb(.)(.„,.)«,...., (39) 

shb(A;)(x, Xo)chb(A;)(?/o, y)Nxydx dy , 
shb(.)(..x,)shb(*)(,..,)Q....d,. 

A similar calculation shows that e^Q^oyoC'^ is a linear combination 

of 

chh{k){x,XQ)chh{k){yo,y)Q^ydx dy , (40) 
shh{k){x,Xo)chh{k){yo,y)N^ydxdy , 
shb(/c)(x, xo)shb(fc)(?/o, y)Q*xydxdy . 

Thus, e^QxoyoQxoyo^~^ is a linear combination of the nine possible 
terms obtained by combining the above. 

Now we list all terms in e^Ve~^Q. Terms in e^Ve~^ ending in 
Qxy are automatically discarded because they contribute nothing when 
applied to Q. The remaining six terms are listed below. 

chb(A;) (xi, xo)chb(/i;)(?/o, yi)shh{k){x2, xo)chh{k){yo, 1/2) 
v{xo - yo)Ql^y^N^2y2^dxi dyi dx2 dy2 dxo dyo , (41) 

chb(A;)(xi, Xo)chb(A;) (yo, yi)shh{k) (xa, Xo)shb(A;) (yo, t/a) 
v{xq - yo)Ql^y^Ql^y^Qdxi dyi dx2 dy2 dxo dyo , (42) 
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shb(A;)(xi, Xo)chb(A;)(yo, yi)shh{k) (xa, Xo)chb(A;)(t/o, I/2) 
v{xo - yo)N^^y^N^^y^Qdxi dyi 6x2 dy2 dxo dyo , (43) 

shb(A;)(xi, xo)chb(fc)(?/o, yi)shh{k){x2, xo)shb(/c)(?/o, 2/2) (44) 
^^(a^o - yo)N^iyiQ*x^y^^dxi dyi dx2 dy2 dxo dyo , 

shb(/c)(xi, Xo)shb(A;)(yo, Z/i)shb(A;)(x2, Xo)chb(A;)(?/o, Z/2) 
^(a^o - yo)QxiyiN^2y2^dxi dyi dx2 dy2 dxo dyo , (45) 



shb(A;)(xi, a;o)shb(A;) {yo, yi)shh{k) {x2, Xo)shb{k){yo, 2/2) 
vixo - yo)QxiyiQl^y^^dxi dyi dx2 dy2 dxo dyo . (46) 

To compute the above six terms, recall f l2^ through ( 1271) as well as 
([T]). In general, N^y^l = l/26{x — y)Q, while / f{x,y)Q*ydxdyfl = 
(0, 0, f{x, y), 0, ■ ■ ■ ) up to symmetrization and normalization. 

The resulting contributions (neglecting symmetrization and normal- 
ization) follow. 

From fHTD: 



Hxuyi)= (47) 
chb(.,(...„)chb(.)fe„,,0*b(^fe,.o)chb(.,(.„...)..(.„ - .„) 
X dx2 dxo dyo . 
From (|42]): 

ipixi,yi,X2,y2) = (48) 

chb(.)(.,,..,chb(.,(.„..OBhb(.,(.„.„,*b(.)feo,.>(.„ - .„) 

X dxo dyo . 
From dlH]): 

i; = (49) 
shb(.)(..,.„)chb(.)fe„,.Oshb(.)(...o)chb(.)fe„,.,).(.„-,„) 

X dxi dx2 dxo dyo . 
From (l44l) . with the and Q* reversed, we get 
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'4'{x2,y2)= (50) 

■shb(.)(x„.o)chb(.)(.o,.Oshb(.)(...o)shb(.)(.o,..).(.o-.o) 
X dxi dxQ dyo , 
as well as the contribution from [A^, Q*], i.e. 

i'(yi,y2)= (51) 

■.hb(.)(...„)chbWfe„,.0*b(.)(..,.„)=hb(i)feo,..).(xo-.„) 

X dxi dxo dyo . 

The contribution of (H5l) is zero, and, finally, the contribution of (146|) . 
using consists of four numbers, which can be represented by the 
two formulas 



= y shh{k){xi, xo)shh{k){yo, Xi)shh{k){x2, xo)shh{k){y(), X2)v{xo — yo) 

(52) 

X dxi dx2 dxo dyo 

and 

ip = \shh{k)\^{xi,Xo)\shh{k)\^{yo,yi)v{xo-yo)dxidyidxodyo . (53) 



We can now state the following proposition. 

Proposition 7.1. The state e^Ve~^Q has entries on the zeroth, sec- 
ond and fourth slot of a Fock space vector of the form given above. In 
addition, if 

II (^i^ - - sh{k)\\Li[o,T]L^dxdy) < Ci, 

,d \ .. 

i-Q^ — + Ay 1 p\\L^[o,T]L^{dxdy) < C2 

and f (x) = or v{x) = xi^)-^^, then 

We^Ve-'^VLf^ dt<C , 



where C only depends on Ci and C2. 
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Proof. This follows by writing ch(A;) = 6{x — y) + p and applying 
Cauchy-Schwartz and local smoothing estimates as in the work of Sjolin 
[32] . Vega [33]; see also Constantin and Saut p]. In fact, we need the 
following slight generalization (see Lemma [7.21 below): If 

II (^-^ - ^xi - ^X2 ± ^X3 ■ ■ ■ ± f{t,Xi, ■ ■ ■Xn)\\L^O,T]L^{dtdx) < C 



with initial conditions 0, then 

f{t, Xi,X2, ■ ■ 



\Xi - X2 



\L^[0,T]L^{dxdy) < C . (54) 



We will check a typical term, (1481) . This amounts to proving the fol- 
lowing three terms are in L^. 

(1) 

ippp{t,xi,yi,X2,y2) = 

pit, xi, Xo)p(t, 2/0, 2/i)shb(A;)(t, Xs, Xo)shb(/c) (t, yo, y2)v{xQ - yo) dxo dyo 



We use Cauchy-Schwartz in xo,yo to get 

\ippp\'^dt dxi dx2 dyi dy2 
< sup j \p{t,Xi,Xo)p{t,yQ,yi)\'^dxidxo dyi dyo 
X / \shb{k){t,X2,xo)shh{k){t,yo,y2)v{xo — yo)\'^dtdx2dxody2dyo< 



JO J 

The first term is estimated by energy, and the second one is 
an application of flMl) with / = shb(A;)shb(A;). Notice that, 
because of the absolute value, we can choose either sh(fc) or 
sh(A;) to insure that the Laplacians in xq, yo have the same 
signs. 

(2) 

il^p5{t,xi,yi,X2,y2) = 

p{t, Xi, Xo)shh{k){t, X2, Xo)shh{k){t,yi,y2)v{xo - yi) dxo . 
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Here, we use Cauchy-Schwartz in xq to estimate, in a similar 
fashion, 

T 

2, 



{ippsl dt dxi dx2 dyi dy2 
< sup J \p(t,Xi,Xo)\^dxidxo 

X / \shh{k){t,X2,xo)shb{k){t,yi,y2)v{xo — yi)\'^dtdx2dxody2dyo<C. 



(3) 

^5s{xi,yi,X2,y2) = shb{k){t, X2, xi)shb{k){t,yi,y2)v{xi - yi) , 

which is just a direct apphcation of (!54l) . 
All other terms are similar. □ 

□ 

We have to sketch the proof of the local smoothing estimate that we 
used above. 

Lemma 7.2. /// : M^n+i ^ £ satisfies 
with initial conditions /(O, ■ ■ ■ ) = 0, then 

f{t, Xi,X2, ■ ■ 



\Xi - X2\ 



L2[0,T]L2(dx) < C . 



Proof. We follow the general outline of Sjolin, [32]. Using Duhamel's 
principle, it suffices to assume that 

d \ 
i— - A,, - A,, ±A,3---±A,J /(t,xi,---x„) = (55) 

with initial conditions /(O, ■ ■ ■ ) = /o G L^. Furthermore, after the 
change of variables Xi ^^^^ i ^2 suffices to prove that 

/()f:,Xi,X2, ■ 



|L2[o,T]L2(dx) 



< C 



Fil 

where / satisfies the same equation fl55l) . Changing notation, denote 
X = {x2,X3, ■ ■ ■ ) and let < ^ >^ be the relevant expression ±|^2p ± 
\^3?.... Write 

f{t,x„x) = I e^*(l«^l'+<«>')e^^'-«^+^^-«/o(6,0 d^id^ ■ 
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Thus, we obtain 
|/(t,xi,x)p 



dtdxidx 



2 

X dtdxdxi 




\Xi 



J \^i-Vi\ 
< j \?o{ii.O?dx,di , 

because one can easily check that 

sup / 5(161' - \vi?)Tr^ — -,dTii < C . 



Thus, the kernel 5(|'CiP~l^iP) is bounded from L'^{drji) to L'^{dC,i). 

1 vi ^ 

□ 

8. Error terms e^[A,V]e~^ 

We proceed to check the operator e^[A,V]e^^ . The calculations 
of this section are similar to those of the preceding section with the 
notable exception of flMl) - flMl) . Recall the calculations of Lemma [3.11 
and write 

e^[A, V]e-^ = j v{x- y) {4>{y)e^ale-^e^a^aye-^ (56) 

+ (f){y)e^ala*ye~^e^a^e'^^ dxdy . (57) 

Now fix xq. We start with the term fl56l) . According to Theorem I4.H 
we have 



^ ^*xo^ = J \^sh{k){x,Xo)ax + ch.{k){x,Xo)alj dx 
while ax^ay^e'^ has been computed in ( l40l) . The relevant terms are 
'shb(.)(x.x„)chb(»:)(,„,y)iV.,<ixd, and 

shb(A;)(a;, xo)shb(/c)(?/o, y)Qlydxdy . 
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Combining these two terms, there are three non-zero terms (which will 
act on Q): 

(1) 



v{xo - i/o)0(l/o)shb(A;)(xi,xo)shb(A;)(x2,xo)shb(A;)(?/o,l/2)axi<5^2j,2f^ 

(58) 

X dxidx2 dy2 dxo dyo . 
This term contributes terms of the form 

i^{t,y2) = j v{xq - yo)(j){t, yo){shh{k){t, xi, xo)fshh{k){t, yo, y2)dxi dxo dyo 

(59) 

as well as the term 

i'{t,X2) = J v{xo - yo)'(p{t,yo)shh{k){t,Xi,Xo)shh{k){t,X2,Xo)shh{k){t,yo,Xi) 

(60) 

X dxi dxo dyo , 

which we know how to estimate. The second contribution is 

(2) 

v{xo - yo)4>{yo)chh{k){xi, Xo)shb(A;)(x2, Xo)chh{k){yo 

(61) 

X dxidx2 dy2 dxo dyo . 

Commuting a*^ with a^^-, we find that fl^T]) contributes 

^(^,2/2) = j v{xo - ?/o)0(t,2/o)chb(/c)(t,Xi,a;o)shb(A;)(t,a;i,Xo)chb(/c)(t,?/o,l/2) 

(62) 

X dxi dxo dyo ■ 

We expand chb(A;)(t, Xi, Xo) = (5(xi — Xq) +p(/i;)(t, Xi — Xq). The 
contributions of p are similar to previous terms, but 5{xi — xo) 
presents a new type of term, which will be addressed in Lemma 
18.21 These contributions are 



'ip5p{t,y2) = J v{xi-yo)(j){t,yo)shb{k){t,xi,xi)p{k){t,yo,y2) (63) 
dxi dyo 

and 

^55(^,2/2) = 4>{t,y2) / v{xi - y2)shh{k){t, xi, xi)dxi . (64) 
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The last contribution of (1561) is 



(3) 

vixo - ?/o)0(l/o)chb(A;)(xi,Xo)shb(A;)(x2,Xo)shb(A;)(?/o, 1/2X1 Q;^2y2^ 
X dxidx2 dy2 dxo dyo ~ 4'{xi,X2, 1/2) 

where 

^p{t,Xl,X2,y2) 

vixo - yo)(f){t, yo)chh{k){t, xi, xo)shb(/c) (t, X2, xo)shb(fc) (t, yo, y2)dxQ dyo 



modulo normalization and symmetrization. This term, as well 
as all the terms in fl57|) . are similar to previous ones and are 
omitted. 

We can now state the following proposition. 

Proposition 8.1. The state e^[A, V]e~^Q has entries in the first and 
third slot of a Fock space vector of the form given above. In addition, 

II - - Ay^ sh{k)\\Li[o,T]L^dxdy) < Ci, 

II " + 1 P\\L^[0,T]L^(dxdy) < Ci 

and 

\\shb{k){t, X, x)\\L2(^[o,T]L^dx)) < C3 , (65) 
and v{x) = for x cl cut-off function, then 

r \\e''[A,V]e-''nf^<C , 
Jo 

where C only depends on Ci, C2, C3. 

Proof. The proof is similar to that of Proposition 17.11 the only excep- 
tion being the terms (l63l) . (l64l) . It is only for the purpose of handling 
these terms that the Coulomb potential has to be truncated, since the 
convolution of the Coulomb potential with the function shb(A;)(a;, x) 
does not make sense. If v is truncated to be in L^{dx) , then we estimate 
the convolution in L'^{dx), and take G L°°{dydt). □ 

□ 

To apply this proposition, we need the following lemma. 
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Lemma 8.2. Let u G X5'3+. Then, 

Proof. As it is well known, it suffices to prove the result for u satisfying 
(^i-^ - A^. - Ay^ u{t, x,y) = 

with initial conditions u{0,x,y) = Uo{x,y) G H^. This can be proved 
as a "Morawetz estimate", see [1^, or as a space-time estimate as in 
[19]. Following the second approach, the space-time Fourier transform 
of u (evaluated at 2^ rather than ^ for neatness) is 

u{t, 20 = c / 6{r - - r]\^ - \^ + r]\^)uo{^ -V,^ + v)dv 



-I 





- r]\ 


\'-\^ + v\ 




i\^-v\ 


+ |e + r^|)V2 



where F{C - V,^ + v) = (1^ - ^1 + 1^ + vlY^^M^ - V,^ + v)- By 
Plancherel's theorem, it suffices to show that 

This, in turn, follows from the pointwise estimate (Cauchy-Schwartz 
with measures) 

~ J + + 

X / Sir-\^-r^\'-\^ + rj\')\F{^-r],^ + r])\'dv 



and the remark that 

-dr] < C 





T]\ 






\^-r]\ + \^ + v\ 





□ 



□ 
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9. The trace / d{t,x,x)dx 

This section addresses the control of traces involved in derivations. 
Recall that 

d{t, X, y) = {ish{k)t + sh{k)g^ + gsh{k)) sh{k) 
- {ich{k)t + [g, ch{k)]) ch(fc) 
— sh(fc)mch(A;) — ch(A;)msh(A;) . 

Notice that if ki{x,y) G L'^{dxdy) and k2{x,y) G L'^{dxdy) then 

J \kik2\{x,x)dx < J \ki(x,y)\\k2{y, x)\dy dx 
< ||^i||l2|I^2||l2 . 

Recall from Theorem 16.11 that if v{x) = -j^ or v{x) = xi^)-^\ then 
ish.{k)t+sh.{k)g'^+gsh.{k), ich.{k)t+[g, ch(fc)] and sh(A;) are in L°°([0, l])L'^{dxdy). 
This allows us to control all traces except the contribution of 6{x — y) 
to the second term. But, in fact, we have 

ich{k)t + [(?, ch(fc)] = {ikt — A^k — Ayk) k — k{ikt — A^k — Ayk) + . . . , 
which has bounded trace, uniformly in [0, 1]. 
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